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ABSTRACT 

We present an analysis of the longest timescale microlensing events discovered by the MA- 
CHO Collaboration during a seven year survey of the Galactic bulge. We find six events that 
exhibit very strong microlensing parallax signals due, in part, to accurate photometric data 
from the GMAN and MPS collaborations. The microlensing parallax fit parameters are used in 
a likelihood analysis, which is able to estimate the distance and masses of the lens objects based 
upon a standard model of the Galactic velocity distribution. This analysis indicates that the most 
likely masses of five of the six lenses are > 1 Mq , which suggests that a substantial fraction of 
the Galactic lenses may be massive stellar remnants. This could explain the observed excess of 
long timescale microlensing events. The lenses for events MACHO-96-BLG-5 and MACHO- 
98-BLG-6 are the most massive, with mass estimates of M/Mq = 6^3° and M/Mq = 6^3, 
respectively. The observed upper limits on the absolute brightness of main sequence stars for 
these lenses are < 1 Lq, so both lenses are black hole candidates. The black hole interpretation 
is also favored by a likelihood analysis with a Bayesian prior using a conventional model for 
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the lens mass function. We consider the possibility that the source stars for some of these six 
events may lie in the foreground Galactic disk or in the Sagittarius (SGR) Dwarf Galaxy be- 
hind the bulge, but we find that bulge sources are likely to dominate our microlensing parallax 
event sample. Future HST observations of these events can either confirm the black hole lens 
hypothesis or detect the lens stars and provide a direct measurement of their masses. Future ob- 
servations of similar events by SIM or the Keck or VLTI interferometers (Delplancke, Gorski 
& Richichi 2001) will allow direct measurements of the lens masses for stellar remnant lenses 
as well. 



1. Introduction 

The abundance of old stellar remnants in our Galaxy is largely unknown because they emit little radia- 
tion unless they happen to be accreting material from a companion star, or for neutron stars, if they happen 
to emit pulsar radiation in our direction. Gravitational microlensing surveys (Liebes 1964; Paczyhski 1986; 
Alcock et al. 1993; Aubourg et al. 1993; Udalski et al. 1993; Bond et al. 2001) have the potential to detect 
completely dark stellar remnants, but for most microlensing events, the mass can only be estimated very 
crudely based upon the observed Einstein ring diameter crossing time, t. For an individual microlensing 
event, the mass can only be estimated so crudely that a 7Mq black hole cannot be distinguished from a 
0.5Mq star. However, for some microlensing events, it is possible to measure other parameters besides t 
that allow tighter constraints on the lens mass (Refsdal 1966; Gould 1992; Nemiroff & Wickramasinghe 
1994; Alcock et al. 1995; Bennett et al. 1996; Han & Gould 1997; Alcock et al. 1997c; Afonso et al. 2000; 
Alcock et al. 2001a). For long timescale microlensing events, which are often due to massive lenses, it is 
frequently possible to measure the microlensing parallax effect (Refsdal 1966; Gould 1992; Alcock et al. 
1995) which is an observable deviation in the microlensing light curve due to the orbital motion of the Earth. 
In this paper, we present an analysis of the microlensing events discovered by the MACHO Project which 
give a very strong microlensing parallax signal, and we show that some of these events are best explained as 
microlensing by black holes. 

This paper is organized as follows: in Sections 2 and 3, we discuss the microlensing event data set and 
the long timescale sub-sample. The microlensing parallax fits are presented in Section 4, and in Section 5 we 
present our main analysis to determine the distances and masses of the lenses. This includes a discussion of 
the projected lens velocity distributions, the source star color magnitude diagrams, and a likelihood analysis 

of the distances and masses of the microlenses. In Section 6, we discuss possible follow-up observations 
with high resolution telescopes and interferometers that can directly determine the microlensing parallax 
event lens masses, and we conclude in Section 7. 
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2. The Data Set 

The MACHO Project (Alcock et al. 1993) has monitored ~ 10-20 miUion stars in the Galactic 
bulge for 6-7 months per year during each of the 1993-1999 Galactic bulge seasons. During the last half 
of 1994, real-time microlensing discovery with the MACHO Alert system became possible (Alcock et al. 
1996). (The OGLE collaboration developed this capability the same year (Udalski et al. 1994b).) This 
development allowed much more accurate photometry of the microlensing events, which were discovered 
in progress, from the CTIO 0.9m telescope where the MACHO/GMAN Project was allocated about 1 hour 
every night (Becker 2000). In 1997, the Microlensing Planet Search (MPS) Project (Rhie et al. 1999) began 
microlensing follow-up observations from the Mt. Stromlo 1 .9m telescope. 

The data set used for this analysis consists of the MACHO survey data from the Mt. Stromlo 1.3m 
"Great Melbourne" telescope for all seven years, CTIO 0.9m data of selected alert events from 1995-1999, 
and Mt. Stromlo 1.9m data of alert events from the MPS 1997-1999 data sets. The initial selection of events 
consists of 42 events from 1993 (Alcock et al. 1997b), 252 events discovered by the MACHO Alert system 
from 1995-1999 (available from http : //darkstar . astro . Washington . edu/), and an additional 
27 events discovered during the testing of the alert system, for a grand total of 321 events. There are 200 
additional Galactic bulge events that have been discovered via other analyses that we have not considered 
here. This paper will focus on the six events from this hst which give a strong microlensing parallax signal. 
The coordinates of these events are given in Table 1. A microlensing parallax study of a larger number of 
MACHO Alert events is presented in Becker (2000). 

The MACHO and MPS data were reduced with slightly different versions of the SoDOPHOT photom- 
etry code (Bennett et al. 1993; Alcock et al. 1999). SoDOPHOT is quite similar to the DOPHOT photometry 

code (Schechter, Mateo, & Saha 1993) that it was derived from, but SoDOPHOT photometry generally ex- 
hibits smaller photometric scatter than DOPHOT photometry. This is due, in part, to SoDOPHOT's error 
flags which allow the removal of suspect data points (Alcock et al. 1999), but the scatter in DOPHOT pho- 
tometry is often increased by the user's choice of PSF fitting parameters. Contrary to expectations, allowing 
the PSF fit box size to scale with the seeing generally causes increased photometric scatter (Bennett et al. 
1993). The photometric errors reported by SoDOPHOT are modified by adding of 1 .4% and 1 .0% in quadra- 
ture to the MACHO and MPS data, respectively, to account for normalization and flat fielding errors. The 
CTIO data were reduced with the ALLFRAME package (Stetson 1994), with the error estimates multiphed 
by a factor of 1.5 to account for systematic errors. 

Table 2 gives the number of observations in each pass band for each event, and the data used for this pa- 
per are presented in Table 3. The complete set of macho survey data is available at http : / / wwwmacho . mcmaster . ca/ 
and http : / /wwwmacho . anu . edu . au. Only the MACHO survey data has been calibrated and trans- 
formed to standard pass bands (Alcock et al. 1999). The other data is given in instrumental magnitudes 
which have only been caUbrated relative to other measurements with the same telescope and passband. 
The transformation between raw MACHO magnitudes given in Table 3 (-Bmacho and -Rmacho) and the 
Kron-Cousins V and R system is given by: 



V = 



2.699 -I- a -Bmacho + & -Rmacho 
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R = 2.412 + cSmacho + ^-Rmacho , 



(1) 



where the coefficients, a, b, c, and d are slightly different for each event as shown in Table 4. 



3. Long Timescale Events 



The timescale of a gravitational microlensing event is described by the Einstein diameter crossing time, 
t, which depends on the lens mass (M), distance (-D^), and transverse velocity {v±). It is given by 



where Dg refers to the distance to the source (typically 8kpc for a bulge source), and Re is the radius 
of the Einstein Ring. Eq. (2) indicates that long t events can be caused by large M, small v± or both. 
Fig. 1 shows the long timescale tail of the t distribution for our sample of 321 Galactic bulge microlensing 
events. In their analysis of the timescale distribution of the 1993 MACHO and OGLE bulge data sets, Han 
& Gould (1996) noted a surprisingly large fraction of the events with t > 140 days: 4/51 or 8%. Such a 
large fraction of long timescale events would be expected less than 2% of the time with any of the stellar 
mass functions that they considered. With our data set of 321 events, we find 28, or 9%, with t > 140 days. 
The MACHO alert system is likely to be somewhat less sensitive to long timescale events than the 1993 
analysis because the alert trigger is based upon the single most significant observation, so we would expect 
a slightly smaller fraction of long timescale events, but the fraction reported here is somewhat higher. The 
formal Poisson probabiUty of 28/321 long events when 2% or less are expected is < 10~^°, so the excess 
of long timescale events over the Han & Gould (1996) models is highly significant. This disagreement may 
be due to a population of massive stellar remnants, including black holes, that was not included in the Han 
& Gould (1996) models, but there are other possibilities as well. Other explanations include a set of source- 
lens systems that have a low relative velocity from our vantage point in the Galactic disk or a more distant 
population of source stars. 

The microlensing parallax effect refers to the effect of the orbital motion of the Earth on the observed 
microlensing light curve. The photometric variation for most microlensing events lasts only a month or 
two. For these events, the change in the Earth's velocity vector during the event is too small to generate a 
detectable deviation from the symmetric light curve, which is predicted for a constant velocity between the 
lens and the Earth-source star fine of sight. For long timescale events, however, it is possible to see the effect 
of the Earth's motion in the microlensing Ught curve, and this is called the microlensing parallax effect. 




(2) 
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4. Microlensing Parallax Fits 

The magnification for a normal microlensing event with no detectable microlensing parallax is given 

by 

A{t) = """ll^ ; u{t) = ,fi[^\^^i^^^2 ^ (3) 

where to is the time of closest approach between the angular positions of the source and lens, and uq = b/RE 
where b is the distance of the closest approach of the lens to the observer-source line. Eq. (3) can be 
generalized to the microlensing parallax case (Alcock et al. 1995) by assuming the perspective of an observer 
located at the Sun. We can then replace the expression for u{t) with 

+2a sm[Q{t — tc)] [uj{t — to) sin + uq cos 6] 

-l-a^sin^/? cos^[i7(t - Q] +2asin/3 cos[ri(t - tc)] [uj{t - to) cos 9 - uq sin 9] (4) 

where A and (3 are the ecliptic longitude and latitude, respectively, 6 is the angle between v± and the North 
ecliptic axis, u = 2/t, and tc is the time when the Earth is closest to the Sun-source line. The parameters a 
and Q, are given by 

a = '^^ilM. (1 _ ecos[Oo(t - t^)]) , (5) 

and 

n{t - tc) = ^o{t - tc) + 2esin[J^o(t - tp)] , (6) 

where tp is the time of perihelion, = 27r yr^^, e = 0.017 is the Earth's orbital eccentricity, and v is the 
lens star's transverse speed projected to the Solar position, given by 

V = v^Ds/{Ds - De) , (7) 



The 28 events shown in Fig. 1 have been fit with the microlensing parallax model described by eqs. (3)- 
(5) which has 5 independent parameters: to, uq, t, v, and 9. In the crowded fields that are searched for 
microlensing, it is also necessary to include two parameters for each independent photometric pass band (or 
telescope) to describe the flux of the source star and the total flux of any unlensed stars that are not resolved 
from the lensed source. Thus, a microlensing parallax fit to the dual-color MACHO data alone will have 9 
fit parameters, and a fit that includes the CTIO and MPS follow-up data will have 13 fit parameters. 

The microlensing parallax fits were performed with the MINUIT routine from the CERN Library, and 
the results for the 6 events that we discuss in this paper are summarized in Table 5. The best fit light curves 
and data are shown in Figs. 2-7. The significance of the microlensing parallax signal is represented by the 
parameter Ax^ shown in Table 5 which is the difference between the fit for a standard microlensing fit 
with no parallax (i.e. v = oo) and the best fit presented here. All 28 events with standard microlensing 
fits (including blending) which indicated t^td > 140 days where fit with a microlensing parallax model 
as well, and the 10 events with a microlensing parallax detection with a significance of Ax^ > 50 are 
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indicated with color in Fig. 1. The four events with 50 < Ax^ < 200 are MACHO-lOl-B, MACHO- 

95-BLG-27, MACHO-98-BLG-1, and MACHO-99-BLG-22, and the six strongest events with Ax^ > 200 
are MACHO-104-C, MACHO-96-BLG-5, MACHO-96-BLG-12, MACHO-98-BLG-6, MACHO-99-BLG- 
1, and MACHO-99-BLG-8. These 6 events are the primary focus of this paper. 

Note that most of the events with t > 200 days and all the events with t > 300 days have a significant 
parallax signal. Microlensing parallax is more easily detected in such long events because the Earth's veloc- 
ity changes significantly during the event, and because long events are likely to have low v± values. There 
are also a number of events with much shorter timescales that appear to have microlensing parallax signals 
significant at the Ax^ > 50 level, but many of these have rather implausible parameters. This is likely to be 
due to the fact that other effects besides microlensing parallax can perturb the microlensing light curves in 
ways that can mimic the parallax effect. Examples of this include binary microlensing (see the discussion 
of MACHO-98-BLG-14 in Alcock et al. 2000a), and the reverse of the parallax effect, the orbital motion 
of a binary source star (Derue et al. 1999; Alcock et al. 2001a; Griest & Hu 1993; Han & Gould 1997), 
sometimes called the "xallarap" effect. The xallarap effect can be particularly difficult to distinguish from 
microlensing parallax because a xallarap light curve can be identical to a parallax light curve if the period, 
inclination, eccentricity, and phase mimic that of the Earth. In practical terms, this is a difficulty only when 
the xallarap or parallax signal-to-noise is weak so that the fit parameters are poorly determined. 

In order to avoid contamination of our microlensing parallax sample with non-parallax microlensing 
events, we have set a higher threshold for the events that we study in detail in this paper: Ax^ > 200. The 
6 events that pass this threshold are listed in Tables 1 and 5. One of these events, MACHO-104-C, was 
the first microlensing parallax event ever discovered (Alcock et al. 1995), and the other five events were 
discovered by the MACHO Alert system. Because of this, they had the benefit of follow-up observations 
by the MACHO/GMAN Collaboration on the CTIO 0.9m telescope or by the MPS Collaboration on the 
Mt. Stromlo 1.9m telescope. Four of these five events would have passed the Ax^ > 200 cut without the 
follow-up data, but event MACHO-98-BLG-6 only passes the cut because of the MPS follow-up data. This 
is probably due to a CCD failure that prevented the imaging of this event in the MACHO-Red band during 
most of the 1998 bulge season. 

We have also compared our microlensing parallax fits to binary lens fits for each of these events. 
The parallax fits are preferred in every case with x^ improvements of 70.3, 81.6, 1625.8, 227.2, 1957.2, 
and I60I.3 for events MACHO-I04-C, MACHO-96-BLG-5, MACHO-96-BLG-I2, MACHO-98-BLG-6, 
MACHO-99-BLG-I, and MACHO-99-BLG-8, respectively. 
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4.1. HST Observations of MACHO-96-BLG-5 

Event MACHO-96-BLG-5 is both the longest event in our sample^ and the event with the faintest 
source star. In fact, the microlensing parallax fit does not constrain the source star brightness very well. 
This is due to the faintness of the source, and due to a potential systematic error. The MACHO camera 
had a CCD upgrade in early 1999 which put a new CCD in the location that views MACHO-96-BLG-5 in 
the MACHO-Red passband. The new CCD probably had the effect of shifting the effective bandpass to a 
slightly different central wavelength, and so a slight systematic shift in the photometry of all the stars might 
be expected to occur with this upgrade. Because the MACHO-96-BLG-5 source is strongly blended with 
unlensed neighbors, the effect of this slight shift on the microlensing fit parameters can be relatively large 
because the fitting routine tries to explain all flux variation as resulting from microlensing. The best fit, 
with this suspect data removed, indicates that only about 12 it 3% of the flux associated with the "star" seen 
in our ground-based images has been microlensed, which would imply that the remaining 88% of the flux 
must come from unlensed neighboring stars which are within ~ 1.5" of the lensed source. Fortunately, we 
have a set of images from the Hubble Space Telescope's WFPC2 Camera that can be used to constrain the 
brightness of the source star more accurately than the fit does. 

We had one orbit of HST data taken in the V and I (F555W & F814W) passbands of the WFPC2 
Camera through Director's Discretionary Proposal # 8490, and this can be used to identify the microlensed 
source star. The first step in this identification process is to determine the centroid of the star that was 
lensed. This can be accomplished by subtracting two images which have substantially different microlensing 
magnifications (Alcock et al. 2001b). Since it is only the lensed source star that will appear to vary in 
brightness, this procedure will yield a point source centered on the location of the lens and source. Of 
course, the subtraction procedure must take into account the differences in the observing conditions of the 
two frames, including differences in seeing, pointing, sky brightness and air mass. We have accompUshed 
this with the use of the DIFIMPHOT package of Tomaney & Crotts (1996). 

A set of 18 of our best CTIO images were selected to use for this source location task because the CTIO 
images generally have better seeing than the MACHO images, and because the highest magnification of the 
source was only observed from CTIO. These 18 images were combined to construct a master reference 
image which was then subtracted from each individual frame to construct a set of 18 difference images. 
The difference frames which had a negative flux at the location of our target were inverted, and then all the 
difference images were combined to make the master difference image shown in Figure 8. The centroid of 
the excess flux in this master difference image can be determined to better than 0.01!' 

In order to identify the lensed source star on the HST images, we must find the correct coordinate 
transformation to match the ground and HST frames, but this is compUcated by the fact that most of the 
"stars" in the ground-based images actually consist of flux from several different stars that are blended 
together in the ground-based frames. We have dealt with this in two different ways: first, we used the HST 

'The analysis of the MACHO and MPS data for MACHO-99-BLG-22 gives a best fit ? = 700 days, but a combined analysis 
with the OGLE data yields a fit that is similar to the OGLE result (Mao et al. 2001) and gives ? = 1100 days. This is about 10% 
longer than our result for MACHO-96-BLG-5. 
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images to select a list of stars that were much brighter than their near neighbors, so that their positions should 
not be greatly affected by blending in the ground-based images. Then, we convolved the HST data with a 
1.2" FWHM Gaussian PSF to simulate the resolution of the ground-based CTIO data. We then analyzed the 
convolved HST image with the same data reduction software used for the ground-based data. This gave an 
additional star hst from the HST image. Two independent coordinate transformations between the ground- 
based and HST data were obtained by matching these two stars hst to the star hst for the ground-based 
data. 

The HST images were dithered, and we combined them with the Drizzle routine (Fruchter & Hook 
1997) prior to the comparison with the ground-based data. The CTIO R-band data were compared to the 

HST 1 and V band images as well as sum of the 1 and V band images. Coordinate transformation were 
determined to match the CTIO image coordinates to each of these HST images using the bright, isolated 
stars in the HST images and with the HST images convolved to ground-based seeing. This resulted in a 
total of 6 different comparisons between the location of the lensed star in the CTIO image and the HST 
images. All these comparisons yielded the same lens star location on the HST frames to better than 0.02", 
and this location coincides with the centroid of the star indicated in Fig. 8. This star was examined carefully 
in both the V and 1 images to determine if it could be a blend of more than one star. Model and DAOPHOT 
generated PSFs were subtracted at the centroid location of the lensed source, but no hint of any additional 
star was found. This star is very likely to be the source star for the MACHO-96-BLG-5 microlensing event. 

The next step in the comparison of the HST and ground-based data is to determine what fraction of the 
flux of the object identified as a star in the ground-based frames is contributed by the source star identified 
in the HST images. This task is comphcated by the fact that there is no close correspondence between the 
passbands of the ground based images and those used for the HST data. (This is due to the limitations 
imposed upon an HST Director's Discretionary time proposal. We requested prompt images in V and I to 
confirm the photometric variation implied by the microlensing parallax model, but prompt imaging in R 
could not be justified. Imaging in R was obtained in a subsequent GO program, and the analysis of these 
data will appear in a future publication.) Presumably, some combination of the V and I band images would 
provide a good representation of the R band ground based image. 

The determination of the lensed flux fraction was made as follows: Photometry of the V-i-I combined 
HST frame was obtained using the IRAF implementation of DAOPHOT (Stetson 1987) and also using 
SoDOPHOT. Both of these packages were also used to reduce the HST images which had been convolved 
to mimic ground based seeing. The total stellar flux of isolated, bright stars was not conserved in these 
convolved images, so we found it necessary to renormalize the stellar flux in the convolved images to the 
ratio found for these isolated bright stars. This comparison yielded a flux fraction of 36% for the lensed 
component of the stellar blend identified as a single star in the ground-based images. We also followed 
this same procedure for the separate I and V images, and the results for the lensed flux fraction were quite 
similar as might be expected from the fit results shown in Table 5, which indicate no color dependence for 
the blending fit parameter. This is likely to be due to the fact that the stars contributing the blended hght and 
the lensed source star are all main sequence stars of similar color which are just below the bulge turn-off. 
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We must also make a correction for the fact that the source star was still being magnified by the lens 
when the HST frames were taken. Because the event timescale depends upon the amount of blending that 
we determine from the HST analysis, it requires an iteration or two to find a fit that predicts the observed 
brightness of the lensed star in the HST frames. The best fit result is that the lensed source provides 33% 
of the total flux of the blended object that would be seen in the ground based frames in the absence of any 
microlensing magnification. At the time of the HST images, the lensing magnification was 1.063 according 
to this fit. 

Finally, we should mention the possibility that the star identified with the lensed source centroid is not, 
in fact, a single star. The HST images reveal no evidence of a chance superposition of unrelated stars, so 
this is unlikely. However, it could be that the superposition is not due to chance. Suppose, for example, that 
the star we've identified as the MACHO-96-BLG-5 source is actually the superimposed images of the lens 
and source. While this is a logical possibility, we will show below that there is no plausible scenario for this 
to occur because the implied lens mass cannot be made compatible with the observed brightness of the lens 
plus source. 



5. Lens Mass and Distance Estimates 

The measurement of the projected speed of the lens, v, allows us to relate the lens mass to the lens and 
source distances 

16G DiDs 16G ,tD, ' ^ ' 

It is often assumed that the distance to the source, Dg, is aheady known, at least approximately, so this 
relation can be considered to give the lens mass as a function of distance. Given the lens distance, one can 
also work out the lens velocity with respect to the line-of-sight to the source, vj,. But, for some distances, 
the imphed v±_ value can be unreasonably small or large. Thus, with some knowledge of the Galactic 
velocity distribution, we can work out an estimate for the distance and mass of the lens. This has been done 
for the MACHO-104-C event using a hkelihood method in Alcock et al. (1995). This analysis assumes 
that the source star resides in the Galactic bulge, which is true for the vast majority of microlensing events 
seen towards the Galactic bulge. The result of similar analyses for the events presented in this paper are 
summarized in Table 6 and in Figs. 11-13. However, the microlensing parallax events are selected from a 
sample of unusually long microlensing events, so it may be that their source star locations are atypical as 
well. With the data currently available to us, we have two ways to investigate the location of the source 
stars for our microlensing parallax events. The first is to make use of the direction of projected velocity as 
determined by the microlensing parallax fit, and the second is to examine the location of the source star in a 
color-magnitude diagram of nearby stars. Another, perhaps more effective, discriminant between different 
source populations is radial velocity measurements. Radial velocities for some of the source stars have been 
measured by Cook et al. (2002), and they have provided us with some preliminary results. 
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5.1. Source Star Locations 

The line-of-sight toward a Galactic bulge microlensing event passes through the Galactic disk, the 
bulge, and through the Sagittarius (SGR) Dwarf Galaxy behind the bulge. So, all of these are possible loca- 
tions for the source stars. The variation in the source population/location can affect the inferred properties 
of the lens in several different ways: 

1. Microlensing rate: The microlensing rate per source star is very much lower for foreground Galactic 
disk stars and very much higher for SGR Dwarf stars than for Galactic bulge stars. Thus, foreground 
disk stars and SGR Dwarf stars will be under-represented and over-represented, respectively, in sam- 
ples of microlensed stars when compared to stars in the Galactic bulge. 

2. Microlensing parallax detectabiUty: some source star populations such as the foreground Galactic disk 
and the Sagittarius Dwarf Galaxy give rise to a larger fraction of events with microlensing parallax 
parameters that can be measured. 

3. Source distance: A source at a greater distance than the nominal Galactic bulge distance will usually 
imply a lower lens mass since M is a decreasing function of Dg in eq. (8) (for fixed x). Similarly, a 
smaller Dg implies a larger mass. 

4. Source velocity: From eq. (7), we see that, for a fixed v value, a smaller v± value implies a smaller Dg 
which, in turn, implies a smaller lens mass (for fixed Di). Smaller v± values are expected for lensing 
of foreground disk sources since the source and lens would both share the Galactic rotation velocity 
of the Sun. 

Several authors who have modeled microlensing parallax events (Mao 1999; Soszynski et al. 2001; 
Smith, Mao, & Wozniak 2001) have suggested that the source stars must be predominantly in the foreground 
Galactic disk because this makes a small v±_ more likely. A disk source is the only possibihty for the OGLE- 
1999-CAR-l event since this star is located far from the bulge, but for events towards the Galactic bulge there 
are several factors that make a foreground disk source star less likely, including a much lower microlensing 
optical depth and a lower density of source stars. These are discussed in section 5.3, where we find that disk 
sources that are definitely in the foreground of the bulge at L>s < 5 kpc are quite unlikely. 

5.2. Projected Velocity Distributions 

One distinguishing characteristic of microlensing parallax distributions for different source populations 
is the distribution of the projected velocity, v including both the amplitude, v, and the direction 0. We use 
a Galactic model in which the stars around us are moving with a velocity dispersion of about 30 km/sec in 
both directions normal to the fine of sight to the bulge. The Sun rotates at a speed of +16 km/sec faster than 
the kinematic Local Standard of Rest (LSR) and is moving towards Galactic North at 7 km/sec (Dehnen & 
Biimey 1998). The Galactic disk rotates with an approximately flat rotation curve at u ~ 200 km/sec, while 
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the Galactic bulge probably has little rotation (Miimiti 1996) and has a velocity dispersion of 80-100 km/sec 
(Spaenhauer, Jones, & Whitford 1992; Minniti 1996; Zoccali et al. 2001). The Sagittarius Dwarf Galaxy is 
moving at 250 ± 90 km/sec in a direction that is only a few degrees away from Galactic North (Ibata et al. 
1997). 

The different velocity distributions of these source and lens populations lead to different expectations 
for the measured v distributions for events from different source star populations. However, the observed 
V distribution is strongly affected by selection effects since only a small fraction of microlensing events 
have detectable parallax signals. These selection effects can be difficult to precisely quantify because of the 
fact that much of the data taken for these events comes from follow-up programs with observing strategies 
that can be subjective and difficult to model. Therefore, instead of attempting a detailed simulation of the 
actual observing conditions, we investigate the v distribution using a "toy model" of a microlensing survey 
and follow-up program. ((Buchalter & Kamionkowski 1997) also performed simulations of microlensing 
parallax events in a somewhat different context.) We assume a disk velocity dispersion of 30km/s in each 
direction, with a flat rotation curve of 200km/s and a bulge velocity dispersion of 80km/s with no bulge 
rotation, and the density profiles are a standard double-exponential disk and a barred bulge as in Han & 
Gould (1996). We assume that events are observed for 7 months per year by a microlensing survey system 
that makes photometric observations with 5% accuracy every 3 days. Once an event is magnified by at least 
0.5 magnitudes, daily follow-up observations start with an accuracy of 1% for each day. This simulated data 
are then fit with a standard, no-parallax microlensing model, and the Ax^ is determined. (Since we have not 
added noise to the light curves, the fit = when there is no microlensing parallax signal.) Events with 

^ 200 are considered microlensing parallax detections, and the v values for these simulated detected 
events are shown in Figure 9. This figure uses Galactic coordinates in which the y-axis is the direction of 
Galactic disk rotation, and the z-axis is Galactic North. 

A striking feature of Figure 9 is that all six of our strong microlensing parallax events have v in the same 
quadrant with positive Vy and negative Vz- This is the region that is preferred for both bulge and SGR source 
stars, but not for foreground disk sources. In our simulations, 65% of the detectable SGR source events, and 
50% of the detectable bulge sources, but only 29% of foreground disk sources lie in this quadrant. 

One selection effect that affects each plot is that events with v roughly parallel to the echptic plane are 
easier to detect than events where v is approximately perpendicular to the ecliptic plane. This effect favors 
the positive -y^-negative Vz and negative -y^-positive Vz quadrants. The reason for this is that the Earth's 
orbital motion only affects u{t) near peak magnification when v is perpendicular to the ecliptic plane, but 
the orbital motion affects u{t) for a longer period of time when it is parallel to v. 

For the bulge sources, there is a preference for positive Vy motion because disk lens stars are passing 
inside of us at a higher angular velocity. If the source stars are rotating with us, as would be the case for 
disk sources in the foreground of the bulge, then the rotation is common to the source, lens and observer, 
and it has no effect. A smaller systematic effect occurs in the disk source case because the Sun is moving 
about 16 km/sec faster than the mean stellar motion around us. Thus, there is a slight enhancement of the 
abundance of negative Vy events. 
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For SGR source stars, signal of the SGR proper motion toward the Galactic North can clearly be seen 

in the strong concentration of events at negative Vz and positive Vy. (Since v is a lens— source velocity, 
the V signal is in the opposite direction of the SGR motion.) For bulge lenses, the 250 km/sec velocity is 
reduced to 90 km/sec by the projection effect, and for the disk lenses that make up the bulk of microlensing 
parallax sample for SGR sources, the typical Vz is —50 km/sec or so. For SGR sources, and disk lenses, the 
combination of SGR proper motion and disk rotation put the majority of v values in the positive Vy-negative 
Vz quadrant where the alignment with the ecUptic plane makes the parallax effect easy to detect. 

5.3. Microlensing Parallax Selection Effects 

For comparison between the different source populations, it is necessary to consider several different 
selection effects. First, the microlensing rate for SGR sources behind our fields is a factor of ~ 6 larger 
than for bulge source stars (Cseresnjes & Alard 2001), and the fraction of SGR events with detectable 
microlensing parallax signals is a factor of 3 larger for sources in SGR than for bulge sources. Thus, it 
would appear that the probability of detecting a microlensing parallax event for a SGR source is a factor of 
~ 20 higher than for a bulge source (assuming that the sources are bright enough for reasonably accurate 
photometry). Of course. Galactic bulge sources are much more numerous, so microlensing parallax events 
with Galactic bulge source stars are likely to be more numerous than events with SGR source stars by an 
amount that is difficult to estimate. We consider this in detail in Section 5.4 when we present the source star 
color-magnitude diagrams. 

It is quite difficult to distinguish Galactic bulge stars from stars in the inner Galactic disk because 
they are at similar distances and their velocity distributions overlap. In fact, this distinction is likely to be 
somewhat artificial because the two components are likely to have merged due to their mutual gravitational 
interactions. Therefore, we will limit our consideration of foreground disk sources to stars with a distance 
< 5kpc. For stars at 5kpc distance at a Galactic latitude of 6 = —3°, the microlensing rate is a factor of 
about ~ 40 lower than for Galactic bulge stars. (The optical depth is only a factor of ~ 20 lower because 
of the longer time scales of disk-disk lensing events.) The physical density of disk stars is about an order of 
magnitude lower than the density of bulge stars, but there is also a volume factor that reduces the number 
of disk stars per unit distance modulus and soUd angle by a factor of 4 at 5kpc. The product of these 
factors yields a net suppression factor of 1/1600 for disk star lensing events for a fixed source star absolute 
magnitude. 

This suppression factor must be multiplied by two enhancement factors. First, our simulations indicate 
that the chances of detecting a microlensing parallax signal are about a factor of 5 larger for disk sources 
than for bulge source stars. This increases the suppression factor to ~ 1/320. There is an additional 
enhancement factor due to the fact that the foreground disk stars are intrinsically fainter and the stellar 
luminosity function rises for fainter stars, but the difference between the disk stars at 5 kpc and the bulge 
stars at 8 kpc is only 1 magnitude. From Holtzman et al. (1998) we see that this factor is at most 10 
if we select a source magnitude such that is 1-2 magnitudes above the bulge main sequence tumoff. For 
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magnitudes that correspond to bulge main sequence stars, it is less than a factor of two. Thus, we expect 

disk stars (with D < 5kpc) to contribute less than 1% of the total number of detectable microlensing 
parallax events, except for source stars that are 1-2 magnitudes above the bulge main sequence turn-off 
where they might account for as many as 3% of the microlensing parallax events with bulge source stars. 

Inner disk stars at D > 5 kpc will be accounted for by allowing their velocities to contribute to the 
assumed bulge velocity distribution. In fact, such inner disk stars are generally not excluded from star 
samples that are used to measure the bulge proper motion (Spaenhauer, Jones, & Whitford 1992; Zoccali et 
al. 2001). We will, therefore, classify all stars in the vicinity of the bulge (5 kpc < Dg < 11 kpc) as bulge 
stars. Instead of trying to distinguish different, but overlapping, populations of source stars, we consider a 
single model including all these stars. 

Stars on the far side of the disk have velocities that make it very unlikely to see the microlensing 
parallax effect, while foreground disk stars are unlikely to be microlensed at all. Therefore, the SGR dwarf 
provides the only "non-bulge" population of potential source stars that we will consider in the remainder of 
this paper. 

5.4. Color Magnitude Diagrams 

Fig. 10 shows color-magnitude diagrams for all the stars within 2 arc minutes around each of our 
microlensing parallax source stars, with the lensed source indicated by a red circle. It is necessary to use 
different color magnitude diagrams for each event because of the large variation in reddening between 
different fields. By plotting only the stars within 2 arc minutes of our targets, we have minimized the 
variation in reddening. 

These CM diagrams indicate that the MACHO-I04-C and MACHO-96-BLG-I2 source stars are lo- 
cated in the bulge red clump region, which means that they are likely to reside in the Galactic bulge. The 
MACHO-99-BLG-8 source star is more luminous than the red clump and is likely to be a bulge giant. Cook 
et al. (2002) find a radial velocity of Vr = 195 it 2km/s which confirms the bulge interpretation for this 
event. The MACHO-96-BLG-5 source star appears to be fainter than virtually all of the other stars in its 
CM diagram. This is a consequence of the extreme crowding of these Galactic bulge fields. The density of 
bright main sequence stars is ~ 2 per square arc second, so main sequence stars are not individually resolved 
in these crowded Galactic fields. Instead, it is groups of unresolved main sequence stars that are identified 
as single stars, and it is these unresolved blends of multiple stars that make up the majority of the fainter 
objects identified as stars in these images. The majority of microlensed source stars in the Galactic bulge 
are blended main sequence stars like the MACHO-96-BLG-5 source, but the microlensing parallax signal is 
easier to detect for brighter source stars. 

The source stars for events MACHO-98-BLG-6 and MACHO-99-BLG-1 appear to be on the bulge sub- 
giant branch of the color magnitude diagram. They have a similar color to bulge red clump stars, but they 
are about 2 magnitudes fainter. This suggests that they could be red clump stars ~ 14 kpc behind the bulge 
in Sagittarius (SGR) Dwarf Galaxy. This is about the only location on the color magnitude diagram were 
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we might expect to see microlensing of SGR source stars, because SGR red clump stars are probably the 
only abundant type of SGR stars that are brighter than the bulge main sequence stars that set the confusion 
limit. This SGR source interpretation appears to gain support from the location of these events in Fig. 9 
which indicates that their parallax velocities are among the ones most consistent with Sagittarius Dwarf 
kinematics. 

A rough estimate of the probability of detecting microlensed SGR source stars can be made by noting 
that SGR Dwarf RR Lyrae stars are about 2.6% as numerous as bulge RR Lyrae in the MACHO fields 
(Alcock et al. 1997a). In a microlensing parallax sample, we should expect SGR source stars to be enhanced 
by a factor of ~ 20, but we must also include both bulge sub-giants and giants in the comparison with the 
SGR red clump giants. This would reduce the fraction of SGR events by a factor two or so. This would 
suggest that we might expect that for every 4 microlensing parallax events with bulge giant or sub-giant 
source we could expect one SGR giant source star event? On the other hand, the ratio of red clump stars to 
RR Lyrae is likely to be higher for SGR stars than for Galactic bulge stars because of the lower metalicity 
of SGR, so we might expect fewer SGR events than this RR Lyrae comparison would suggest. 

These considerations suggest that we should take the SGR source star hypothesis seriously for these 
events. However, Cook et al. (2002) used the Keck HIRES spectrograph to obtain spectra of the source stars 
for these events, and they find radial velocities of Vr = —65 ± 2 km/s and = 64 ± 2 km/s for MACHO- 
98-BLG-6 and MACHO-99-BLG-1, respectively This is not consistent with the SGR radial velocity (Ibata 
et al. 1997) of Vr = 140 ib 10 km/s, and they are about 2a away from the expectation for a disk source star 
(Wielen 1982). Thus, these events are most likely to have bulge sub-giant source stars. 

5.5. Likelihood Distance and Mass Estimates 

Another, somewhat more general, constraint on x and M can be obtained if we make use of our knowl- 
edge of the velocity distributions of the source and lensing objects, since the likelihood of obtaining the 
observed value of v is a strong function of the distance to the lens. Note that this assumes that stellar 
remnant lenses have a velocity and density distribution that is similar to that of observed stellar popula- 
tions. For neutron stars, this might be a questionable assumption because many neutron stars are apparently 
bom with a large "kick" velocity. However, for black holes, the evidence indicates that significant kick 
velocities are rare (Nelemans, Tauris, & van den Heuvel 1999). As an example of such an analysis, let us 
suppose that the disk and bulge velocity dispersions were negligible relative to the Galactic rotation velocity. 
Then, for disk lenses we would obtain the relation v = 200D(/{Ds — D() km/s implying a lens distance 
of = Dgv/ {v + 200 km/s). In reaUty, the random motions of both disk and bulge stars broaden this 
relationship somewhat, but we can still obtain a useful constraint. 

^ A previous estimate of lensing rates for SGR source stars has been made by Cseresnjes & Alard (200 1 ) who find a smaller ratio 
of SGR/bulge source lensing events than our estimate. This is because they do not consider only microlensing parallax events and 
because they count the much more numerous events with bulge main sequence source stars. 
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Given the observed v, we obtain a likelihood function 

L{x; v) oc - a;) pl{x) v{l - xf J fsi^s) - x){vq + v) + xvs) dvs, (9) 

where is the density of lenses at distance x = D^/Ds, and the integral is over combinations of source 
and lens velocities giving the observed v. V5 and = (1 — x)(vo + v) + xys are the 2-D source and 
lens velocity distribution functions (normalized to unity). We assume the same Galactic parameters as in 
our V simulations above: a disk velocity dispersion of 30 km/s in each direction, a flat disk rotation curve 
of 200 km/s, and a bulge velocity dispersion of 80 km/s with no bulge rotation. The density profiles are 
a standard double-exponential disk and a Han & Gould (1996) barred bulge. For all events, the source is 
assumed to reside in the bulge, while the lens may be in the disk or the bulge. But for events MACHO-98- 
BLG-6 and MACHO-99-BLG-1, we also consider the possibiUty of a SGR Dwarf source star with the lens 
in the disk or bulge, although this now appears to be ruled out (Cook et al. 2002). 

The resulting likelihood functions for is shown as the long-dashed curves in Figs. 11-13, and these 
are insensitive to specific parameter choices. These likelihood functions also provide a means for estimating 
the lens masses via the relation (8), which is also plotted in Figs. 11-13. Fig. 14 shows how the mass 
estimates correlate with the best fit event timescale for the six high signal-to-noise microlensing parallax 
events as well as four other events of lower signal-to-noise. (The lower signal-to-noise event with the highest 
mass is MACHO-99-BLG-22/OGLE-1999-BUL-32 which has been presented as a black hole candidate by 
Mao etal. 2001.) 

One common way to interpret likelihood functions is the Bayesian method, in which the lens mass 
(or distance) probability distribution is given by the likelihood function times a prior distribution, which 
represents our prior knowledge of the probability distribution. In our case, the hkeUhood function represents 
all of our knowledge about the lens mass and location, so we select a uniform prior. With a uniform prior, 
the Ukelihood function becomes the probability distribution and we are able to calculate the lens mass 
confidence levels Usted in Table 6. This table also includes lens mass confidence levels for models that 
differ from the preferred model in order to show how the mass estimates depend upon the amount of blending 
(for MACHO-96-BLG-5) and on whether the source star resides in the Galactic bulge or the SGR Dwarf. 
Note that the uncertainty in the mass estimates is smaller for SGR Dwarf sources due to the small velocity 
dispersion of the SGR Dwarf and the smaller range of likely lens distances. 



5.6. Constraints on Main Sequence Lenses 

If we assume that the lens stars are main sequence stars, then we can obtain an additional constraint on 
their distances and masses by comparing the brightness of a main sequence star, of the implied mass, to the 
upper limit on the brightness of the lens star. We have assigned a conservative upper limit on the V-band 
brightness of each lens star based upon the available photometry and microlensing parallax fits hsted in 
Table 5. In the case of MACHO-96-BLG-5, the upper hmit is particularly stringent because it is based upon 
HST observations. Note that if we assign some of the flux of the star identified in the HST images to the 
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lens star instead of the source, the best fit t will increase almost hnearly with the inverse of the source star 
flux. This causes the lens mass estimate to increase as ~ Since stellar luminosity varies as a high power 
of the mass, a main sequence lens will be more strongly ruled out. 

In order to apply these constraints to the likelihood functions for the mass and distances of the lens 
stars, we have multiplied the likelihood function by the Gaussian probability that the lens brightness ex- 
ceeds the upper limit on the brightness of lens star. If a main sequence lens star would be fainter than the 
observed maximum brightness, there is no modification of the hkehhood function. This gives the short- 
dashed hkehhood curves shown in Figs. 11-13. These results are insensitive to our assumed L oc 
mass-luminosity relation. The assumed maximum lens brightnesses are F = 19.88, 20.57, and 16.92 for 
events MACHO-96-BLG-12, MACHO-98-BLG-6, and MACHO-99-BLG-8, respectively These are based 
upon the amount of blending allowed by the fit, and each of these has an assumed 25% uncertainty which 
is also based upon the fit. For MACHO-96-BLG-5, the maximum lens brightness is F = 23.63, with an 
assumed 50% uncertainty. For MACHO- 104-C, and MACHO-99-BLG-8, the best fit has very httle blended 
flux: V = 22.33, and 23.08, respectively. But, in both cases, the uncertainty in the blended flux is five times 
the best fit value. 

The properties of the most likely main sequence lens models are given in Table 7, which is discussed 
in more detail in section 6.2. An important parameter in this table is the predicted lens-source separation 
in June, 2003, when they might plausibly be observed by HST. This can be calculated from the lens-source 
proper motion which is related to the projected velocity by /x = v{Ds — Dg) j (DgDi). 

5.7. Stellar Remnant Lenses and Black Hole Candidates 

The mean mass estimate for the six microlensing parallax events is 2.7M0. Five of the six have best fit 
masses > 1 M©, and two of the events, MACHO-96-BLG-5 and MACHO-98-BLG-6, have best fit masses 
> 3M0. This makes them black hole candidates because the maximum neutron star mass is thought to be ~ 
2M0 (Akmal, Pandharipande, & Ravenhall 1998). The 95% confidence level lower limits on the masses of 
these lenses are 1.64M0 and O.94M0, respectively, while the 90% confidence level lower hmits are 2.3M0 
and 1.9M0. A main sequence star lens at the lower limit mass is strongly excluded in the case of MACHO- 
96-BLG-5 because of the constraint on the lens brightness from HST images. However, a main sequence 
lens with a mass at the 95% confidence hmit is not quite excluded for MACHO-98-BLG-6. The masses 
that have been measured for neutron stars are close to the Chandrasekhar mass, M^s = 1.35 ± O.O4M0 
(Thorsett & Chakrabarty 1999), which is excluded at better than 95% confidence for MACHO-96-BLG-5 
and better than 90% confidence for MACHO-98-BLG-6. Thus, both MACHO-96-BLG-5 and MACHO- 
98-BLG-6 are both black hole candidates, but there is a small chance that MACHO-98-BLG-6 could be a 
neutron star or even a main sequence star. 

In addition to these black hole candidates, three of the remaining four microlensing parallax events 
have best fit masses > IM0. For MACHO-104-C and MACHO-96-BLG-12, main sequence lens are dis- 
favored, but not ruled out. MACHO-99-BLG-8 appears to be blended with a relatively bright source, so 
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a main sequence lens of M ^ IM© is a possibility. As we explain below, with HST imaging it will be 
straightforward to detect the lenses if they are main sequence stars. If HST images fail to detect the lens 
stars, then we can show that the lenses are almost certainly stellar remnants. 

5.8. Likelihood Analysis with a Mass Function Prior 

The Likelihood analysis presented in Section 5.5 attempts to estimate the distance to the lens based 
upon the measured value of the projected velocity, v, and then the lens mass is determined from eq. 8. If 
the lens mass function, dn/dM = 4>{M), is known, then it is possible to use the measured t value to make 
a more accurate estimate of the lens mass as advocated by Agol et al. (2002). The UkeUhood function, 
eq. 9, can be modified by multiplying by 5{t — tm)M^/'^(j){M)dM and integrating over M, where tm is the 
measured value of t. The factor of M^/^ is the contribution of the lens mass to the lensing cross section, 
which is proportional to Re- The integral over 5{t — t„i)dM gives an additional factor of M. Thus, the 
likelihood analysis presented in Section 5.5 is equivalent to assuming a mass function of 4>{M) oc M~^'^. 

A more conventional mass function for the Galactic bulge is a broken power law initial mass function 
(Kroupa 2002) with ^(M) oc M'^-^ for O.O3M0 < M < CSM©, and 4>{M) oc M'^-^^ for O.8M0 < 
M < 1OOM0. However, the stars with M > I.OMq will generally have ended their main sequence lifetimes 
and have become stellar remnants after significant mass loss. Following Fryer & Kalogera (2001), we can 
assume that all stars with an initial mass greater than a particular cutoff mass, Mj > Mbh become black 
holes. We take Mbh = 2OM0 (Fryer 1999; Fryer & Kalogera 2001). Similarly, we assume that all stars with 
8M0 < Mi < Mbh become neutron stars, and all stars with l.OM© < Mi < 8Mq become white dwarfs. 
The mass functions of the stellar remnants are assumed to be Gaussians with mean masses of O.6M0, 
I.35M0, and 8M0 for white dwarfs, neutron stars, and black holes respectively. The Gaussian sigmas 
are O.15M0, O.O4M0, and 2.5M0, respectively. These are consistent with the measured mass functions 
(Bergeron, Saffer, & Liebert 1992; Bergeron, Leggett, & Ruiz 2001; Thorsett & Chakrabarty 1999; Bailyn, 
Jain, Coppi, & Orosz 1998), although the difficulty of directly observing old stellar remnants assures that 
the observed samples are incomplete. With this mass function, black holes would account for 3.7% of the 
Galaxy's stellar mass. 

A Bayesian analysis based upon this mass function gives a probability of 93% that the MACHO-96- 
BLG-5 lens is a black hole and a probability of 69% that the MACHO-98-BLG-6 lens is a black hole. The 
probability of at least one black lens is 98%. This analysis may underestimate the black hole probability 
because the assumed mass function caimot account for the large number of long timescale microlensing 
events. An initial IMF that is slightly shallower than the Salpeter slope, (t){M) oc M~'^-^, might be appro- 
priate if most of the stars in the Galaxy were formed in denser or more metal poor regions than is typical for 
present day star forming regions (Figer et al. 1999; Smith & Gallagher 2001). With this mass function and 
with Mbh = 2OM0, black holes would account for 12% of the Galaxy's stellar mass. When we repeat the 
likelihood analysis with this mass function, we find black hole probabilities of 97% for MACHO-96-BLG-5 
and 88% for MACHO-98-BLG-6. The probability of at least one black hole lens with this mass function is 
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99.7%. If we retain the Salpeter IMF slope, and increase Mbh to 4OM0, then the black hole probabilities 
for MACHO-96-BLG-5 and MACHO-98-BLG-6 drop to 82% and 43%, respectively. However, such a mass 
function probably cannot explain the excess of long timescale events. 

We should note that these probabilities are substantially larger than those reported in a similar analysis 
in a preprint by Agol et al. (2002). This was due to a likelihood function calculation error by Agol et al. 
(2002). When this error is corrected, their results are quite similar to those presented here (Agol, private 
communication) . 



6. Follow-up Observations 

The detection of the microlensing parallax effect allows us to make a lens mass estimate that is accurate 
to about a factor of two, and to identify the black hole candidates. However, these estimates are not accurate 
enough to determine the black hole mass function, and they do not allow the unambiguous identification 
of neutron star or white dwarf lenses. However, follow-up observations with higher resolution instruments 
hold the promise of much more precise determinations of the lens masses. 



6.1. Interferometric Follow-up 

The most ambitious of microlensing event follow-up plans involve interferometric instruments such 
as the Keck and VLT interferometers (Delplancke, Gorski & Richichi 2001) and the Space Interferometry 
Mission (SIM) (Boden, Shao, & van Buren 1998). The most spectacular confirmation of a black hole event 
would be to measure the image splitting which is given by 

0sep = 2e£Vl + ^^V4, (10) 

where 6e is the image separation and u is given by eq. (4). For MACHO-96-BLG-5, we have 6e = 9.8 mas 
if the lens is at the distance preferred by the Ukehhood analysis. This compares to the 5 mas diffraction limit 

of an interferometer with a 100 m baseline operating at a wavelength of 2 /im, such as the Keck or VLT 
Interferometers. In fact, these instruments are expected to be able to measure image splittings as small as 
~ 30 //as (Delplancke, Gorski & Richichi 2001). Such measurements would allow a direct measurement of 
the lens mass: 

The most challenging aspect of such measurements is the faintness of source stars such as the MACHO- 
96-BLG-5 source, which is close to the (rather uncertain) magnitude limit of the VLT Interferometer 
(Delplancke, Gorski & Richichi 2001). 



Even if the images cannot be resolved, it may be possible to measure the deflection of the image 



-19- 



centroid (Hog, Novikov, & Polnarev 1995; Miyamoto & Yoshi 1995; Walker 1995) which is given by 

This can be measured by a very accurate astrometry mission such as SIM (Boden, Shao, & van Buren 1998; 
Paczynski 1998; Gould 2000). Once again, however, the MACHO-96-BLG-5 source is a rather faint target 
for SIM, but in this case, the measurement is not so difficult because the ampUtude of the centroid motion is 
very much larger than SIM's sensitivity limit. 

If it should turn out that some of the more massive lenses are located very close to us, then it might 
be possible to directly observe the lensed images with HST. This is a realistic possibility for the MACHO- 
99-BLG-22/OGLE-1999-BUL-32 event (Mao et al. 2001) because its v value is in the opposite quadrant 
from the events studied in this paper. This gives a likelihood function with two peaks: one at a distance of 
~ 500 pc for a lens in the disk and one at a distance of ~ 6kpc for a bulge lens (Bennett et al. 2002). The 
bulge lens solution predicts a mass of a few M0, but the disk lens solution predicts a mass of > IOOM0 and 
a lensed image separation of ~ 0.1'.' 



6.2. Lens Detection and Source Proper Motion 

Another method can be used to make a direct determination of the lens mass for a bright lens star. If the 
lens can be detected and the relative proper motion of the lens with respect to the source is measured, then 
it is also possible to determine the lens mass from the proper motion and microlensing parallax parameters 
with the following formula: 

where /x is the relative lens— source proper motion. This technique has the advantage that the proper motion 
measurements can be made many years after the peak magnification of the microlensing event. The lens- 
source separation can reach the 50-100 mas range within 5-10 years. Table 7 shows the predicted separations 
and lens brightness contrasts for our six strong microlensing parallax events. The columns are (1) the 
MACHO event name, (2) the lens mass with la errors, (3) a Ukely lens mass, MrmMS^ if the lens is on 
the main sequence, (4) the lens distance, -D^_ms for a lens of mass M^mMS' (5) the predicted lens-source 
separation in June, 2003, (6) the apparent V magnitude of the lenses, and (7-11) the predicted contrast 
between the lens and source brightness in the UBVl bands. Positive A-mags. imply that the source is 
brighter than the lens, so lens detection is easiest for events that have small or negative A-mag. values. With 
the exceptions of MACHO-96-BLG-5, which doesn't have a viable main sequence lens model, all of the 
other lens stars should be detectable if they are not stellar remnants. 

When the lens can be detected, it should also be possible to constrain the unlensed brightness of the 
source star, which will reduce the error bars on t. Also, it should be possible to get very accurate measures 
of the relative proper motion, /x, as the lens moves further from the source. Thus, the ultimate Umits on the 
masses of the lenses may come from the uncertainties in the v values, which range from 2 — 10 %. 
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When the lenses are undetectable, it should still be possible to measure the proper motion of the source 
star with HST images separated by ~ 5 years. The proper motion can only be measured with respect to 
the average of other, nearby stars because extra-galactic reference sources are not easily identified in these 
crowded Galactic bulge fields (Spaenhauer, Jones, & Whitford 1992; Zoccali et al. 2001). Proper motion 
measurements of the microlensed source stars would allow us to remove one degree of freedom from our 
UkeUhood analysis and reduce the uncertainty in the impUed lens distances and masses. The proper motion 
distribution of the stars in the same field will also allow us to test the Galactic models that are used for 
the likelihood analysis, and so this should reduce the systematic uncertainties in the lens distance and mass 
estimates. 



7. Discussion and Conclusions 

We have performed microlensing parallax fits on the Galactic bulge events detected by the MACHO 
Collaboration with timescales of i > 140 days, and found six events with highly significant detections of the 
microlensing parallax effect. Our analysis of the velocity distributions expected for parallax microlensing 
events from different source star populations suggests that source stars in the SGR Dwarf Galaxy might 
contribute to the detectable microlensing parallax events, and inspection of the source star color-magnitude 
diagrams indicates that two of our microlensing parallax events have source stars which could be SGR Dwarf 
red clump stars. However, radial velocity measurements (Cook et al. 2002) indicate that they are probably 
bulge sub-giant stars. 

A likelihood analysis has been employed to estimate the distance and masses of the lenses, and this 
indicates an average mass for our six lenses of 2.7 M©. Two of the lenses have masses large enough to 
imply that they are probably massive stellar remnants: The mass estimates for the MACHO-96-BLG-5 and 
MACHO-98-BLG-6 lenses are M/M© = 6+3° and M/M© =6+3, respectively, which impUes that both 
are likely to be black holes. Together with MACHO-99-BLG-22/OGLE-1999-BUL-32 (Mao et al. 2001), 
these are the first black hole candidates that are truly black since we have not seen any radiation from matter 
that is gravitationally bound to the black hole. 

Our likelihood analysis differs from that of Agol et al. (2002) in that we compute the likelihood for 
the measured v value whereas Agol et al. (2002) attempt to compute the likelihood of the measured values 
of t as well as v. However, this requires that we input the mass function of the lenses, and this has never 
been measured for a complete sample of stellar remnants. Thus, the method of Agol et al. (2002) can give 
misleading results if the input mass function is not correct. Nevertheless, the results of such an analysis are 
consistent with the results that we have presented here. (Note that the preprint version of Agol et al. (2002) 
claimed an inconsistency with our results, but this was due to an error in the computation of the likelihood 
function (Agol, private communication).) For the MACHO-99-BLG-22/OGLE-1999-BUL-32 event, the 
method of Agol et al. (2002) does give potentially misleading results, however, because the shape of the 
Likelihood function for this event makes the results quite sensitive to the assumed black hole mass function 
(Bennett et al. 2002), which is, of course, unknown. 
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Similar events detected in the next few years may yield lens masses that are measured much more pre- 
cisely due to follow-up observations from ground-based (Delplancke, Gorski & Richichi 2001) and space- 
based (Gould 2000) interferometers. This will allow an unambiguous determination of the abundance and 
mass function of black hole and neutron star stellar remnants, although it may be difficult to determine if 
~ 2 Mq objects are black holes or neutron stars. At present, there are three black hole microlens candidates 
in the sample of 321 microlensing events that was the starting point for this paper (although MACHO-99- 
BLG-22 is only identified as a strong black hole candidate when OGLE data are included in the analysis 
(Mao et al. 2001; Bennett et al. 2002)). This is about 1% of the events, but far more than 1% of the total con- 
tribution to the microlensing optical depth. This suggests that the fraction of our Galaxy's stellar mass that 
is in the form of black holes may be significantly larger than 1%, which might help to explain the observed 
excess of long timescale microlensing events. However, we have not made an accurate determination of our 
microlensing event detection efficiency for this data set, and the detection efficiency is certainly larger for 
long timescale microlensing events than for short events. It is also possible that one of these three lenses 
may not be a black hole, and so these microlensing results may still be consistent with models which predict 
that of order 1% of the Milky Way's stellar mass should be in the form of black holes (Brown & Bethe 
1994; Fryer & Kalogera 2001; Gould 2000). If all three of these events are truly due to black hole lenses, 
then a black hole mass fraction as high as ~ 10% might be preferred. These results appear to indicate that 
most stellar mass black holes do not reside in the X-ray binary systems where they are most easily observed 
(Bailyn, Jain, Coppi, & Orosz 1998). 
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Fig. 1. — The distribution of event timescales, t, for the 28 events with t > 140 days. The colored bars 
indicate the events with formally significant detections of microlensing parallax, but parallax signal for the 
events indicated in yellow is weak enough that the detection is not considered to be definitive. 
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Fig. 2. — MACHO- 104-C light curves normalized to the unlensed flux of the lensed star. The MACHO red 
and blue data are plotted in magenta and blue, respectively. The black curve is the parallax fit while the cyan 
curve is the best fit standard microlensing lightcurve. An additional 5 years of data showing no photometric 
variation are not shown. 
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Fig. 3. — MACHO-96-BLG-12 lightcurve closeup with lightcurves normalized to the unlensed flux of the 
lensed star. The MACHO red and blue data are plotted in magenta and blue, respectively, and the CTIO data 
are shown in red. The black curve is the parallax fit while the cyan curve is the best fit standard microlensing 
lightcurve. An additional 5 years of data showing no photometric variation are not shown. 
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Fig. 4. — MACHO-98-BLG-6 lightcurve closeup with lightcurves normalized to the unlensed flux of the 
lensed star. The MACHO red and blue data are plotted in magenta and blue, respectively, the CTIO data are 
shown in red, and the MPS data are shown in green. The black curve is the parallax fit while the cyan curve 
is the best fit standard microlensing lightcurve. The gap in the MACHO red data during the day 2280-2650 
interval is due to a CCD failure. An additional year of data showing no photometric variation is not shown. 
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Fig. 5. — MACHO-99-BLG-1 light curves normalized to the unlensed flux of the lensed star. The MACHO 
red and blue data are plotted in magenta and blue, respectively, the CTIO data are shown in red, and the 
MPS data are shown in green. The black curve is the parallax fit while the cyan curve is the best fit standard 
microlensing lightcurve. An additional 4 years of data showing very little photometric variation are not 
shown. 
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Fig. 6. — MACHO-99-BLG-8 light curves normalized to the unlensed flux of the lensed star. The MACHO 
red and blue data are plotted in magenta and blue, respectively, the CTIO data are shown in red, and the 
MPS data are shown in green. The black curve is the parallax fit while the cyan curve is the best fit standard 
microlensing lightcurve. An additional 3 years of data showing very little photometric variation are not 
shown. 
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Fig. 7. — MACHO-96-BLG-5 lightcurves normalized to the unlensed flux of the lensed star. The MACHO 
red and blue data are plotted in magenta and blue, respectively, and the CTIO data are shown in red. The 
black curve is the parallax fit while the cyan curve is the best fit standard microlensing lightcurve. An 
additional 4 years of data showing very little photometric variation are not shown. 



Fig. 8. — The image on the left is the master difference image as described in the text. It has been registered 
to the same coordinate system as the F814W HSTAVFPC2 image shown on the right. The red marks 
show the centroid of the variable flux in the master difference image and the location of this centroid when 
transformed to the coordinate system of the HST data. A single, main sequence bulge star is clearly identified 
as the lensed source star. 
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Fig. 9. — The distribution of v values in Galactic coordinates is shown for simulated microlensing parallax 
events towards the Galactic bulge for three different source star populations: the Galactic bulge, the fore- 
ground Galactic disk, and the Sagittarius Dwarf Galaxy. The large colored dots show the locations of our 
detected microlensing parallax events. The red circular spot is our best black hole candidate, MACHO-96- 
BLG-5, and the green circular disk is the other black hole candidate: MACHO-98-BLG-6. The two green 
spots are the events with source stars that appear to be bulge sub-giants or Sagittarius Dwarf red clump stars. 
(MACHO-99-BLG-1 is the other). The blue squares are the bulge red clump source star events, and the blue 
triangle is MACHO-98-BLG-8 which has a red giant source and is probably also in the bulge. 
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Fig. 10. — Color-Magnitude diagrams from MACHO data are shown for all the detected stars within a 2 arc 
minute circle around each of our microlensing parallax events. The red circles indicate the location of the 
lensed source star, assuming the blending implied by the microlensing parallax fits. 
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Fig. 11. — The mass vs. distance relations (solid curves) for our candidate black hole lenses are shown along 
with the likelihood functions (long dashed curves) computed assuming a standard model for the Galactic 
phase space distribution. The source star is assumed to reside in the bulge for both events. The impUed best 
fit masses are M = etg^M© for the MACHO-96-BLG-5 lens and M = e+^M© for the MACHO-98-BLG- 
6. The 95% confidence level lower limits on the masses are 1.6Mq and O.94M0 respectively. The short 
dashed curves delineate the portion of the likelihood functions that is allowed when the lens is assumed to 
be a main sequence star. The ratio of the area below this portion to the entire area below the hkehhood curve 
gives a probability that a lens is a main sequence star. For MACHO-96-BLG-5, the upper limit on the lens 
brightness is very stringent because of the HST images, and a main sequence lens is ruled out. 
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Fig. 12. — The mass vs. distance relations (solid curves) for our two bulge clump giant source events are 
shown along with the likelihood functions (long dashed curves) computed assuming a standard model for 

the Galactic phase space distribution. The implied best fit masses are M = 1.1+J;^M0 for the MACHO- 
104-C lens and M = for the MACHO-96-BLG-12. The 95% confidence level lower limits 

on the masses are O.35M0 and O.SSM© respectively. The short dashed curves delineate the portion of the 
UkeUhood functions that is allowed when the lens is assumed to be a main sequence star, and they indicate 
that main sequence lenses are disfavored but not ruled out. 
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Fig. 13. — The mass vs. distance relations (solid curves) for the two 1999 microlensing parallax events are 
shown along with the UkeUhood functions (long dashed curves) computed assuming a standard model for the 
Galactic phase space distribution. For both events the source star is assumed to reside in the Galactic bulge. 
The implied best fit masses are M = O.T+q^Mq for the MACHO-99-BLG-1 lens and M = 1.2+ggM0 
for the MACHO-99-BLG-8. The 95% confidence level lower Hmits on the masses are O.14M0 and 0.3Mq 
respectively. The short dashed curves dehneate the portion of the likelihood functions that is allowed when 
the lens is assumed to be a main sequence star, and these indicate that the lens brightness constraints are 
consistent with main sequence lens stars. For MACHO-99-BLG-1, a main sequence lens is disfavored, 
however. 
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Fig. 14. — This plot shows M vs. t for the 10 events with t > 140 days with la error bars for the mass 
estimates. All events with formally significant detections of microlensing parallax are shown, but parallax 
signal for the events indicated in yellow is weak enough that the detection is not considered to be definitive. 
The green open symbols indicate the predicted lens masses for MACHO-98-BLG-6 and MACHO-99-BLG- 
1 if their source stars were in the SGR Dwarf Galaxy, a possibility that appears to be contradicted by their 
spectra. 
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Table 1. Microlensing Parallax Event Coordinates 



Galactic Ecliptic 



Event Name 


MACHO Star ID 


RA(J2000) 


DEC(J2000) 


I 


b 


A 


P 


MACHO- 104-C 


104.20251.50 


18 


03:34.0 


-28:00:19 


2.797 


-2.933 


270.790 


-4.568 


MACHO-96-BLG-5 


104.20906.3973 


18 


05:02.5 


-27:42:17 


3.219 


-3.071 


271.119 


-4.270 


MACHO-96-BLG-12 


104.20382.803 


18 


03:53.2 


-27:57:36 


2.871 


-2.973 


270.861 


-4.524 


MACHO-98-BLG-6 


402.48103.1719 


17 


57:32.8 


-28:42:45 


1.526 


-2.132 


268.762 


-5.267 


MACHO-99-BLG-1 


121.22423.1032 


18 


08:50.0 


-30:31:56 


1.138 


-5.162 


271.917 


-7.106 


MACHO-99-BLG-8 


403.47849.756 


17 


56:25.2 


-29:40:31 


0.569 


-2.401 


269.218 


-6.237 



Table 2. Number of Observations 



Event 


MACHO-Red 


MACHO-Blue 


CTIO 


MPS 


104-C 


534 


308 








96-BLG-5 


558 


1542 


179 





96-BLG-12 


584 


466 


103 





98-BLG-6 


952 


1083 


29 


212 


99-BLG-l 


343 


260 


11 


153 


99-BLG-8 


386 


310 


213 


155 
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Table 3. Photometric Measurements 



Event Name Pass Band time (MJD) Magnitude uncertainty 



MACHO-104-C MACHO-Red 



MACHO-Blue 



MACHO-96-BLG-5 MACHO-Red 



MACHO-Blue 



CTIO 



430.79500 


14.0310 


0.0221 


438.78620 


13.9090 


0.0188 


441.73940 


13.8060 


0.0163 


442.74640 


13.8120 


0.0163 


443.71500 


13.7760 


0.0172 


446.72880 


13.7400 


0.0155 


453.79520 


13.5360 


0.0182 


459.71700 


13.3230 


0.0301 


463.67350 


13.2520 


0.0172 


463.67660 


13.2570 


0.0200 


430.79500 


13.0240 


0.0167 


AA2.1A(A0 


12.7870 


0.0157 


443.71500 


12.7820 


0.0157 


452.74020 


12.5650 


0.0157 


453.79520 


12.5320 


0.0160 


455.75470 


12.4810 


0.0157 


457.78640 


12.4300 


0.0157 


459.71700 


12.3510 


0.0200 


463.67350 


12.2290 


0.0157 


463.67660 


12.2500 


0.0167 


430.79500 


16.2340 


0.1299 


441.73940 


16.3030 


0.0726 


AA2.1AM0 


16.4470 


0.0814 


443.71500 


16.2560 


0.0942 


455.75470 


16.2170 


0.0952 


457.78640 


16.3910 


0.1110 


459.71700 


15.8800 


0.2295 


463.67350 


16.3050 


0.1538 


463.67660 


15.9980 


0.1796 


465.65960 


16.4010 


0.2574 


430.79500 


17.0900 


0.2155 


438.78620 


17.2230 


0.2066 


441.73940 


17.1900 


0.1031 


442.74640 


17.3520 


0.1140 


443.71500 


17.2020 


0.1388 


446.72880 


17.4260 


0.0580 


452.74020 


17.1390 


0.1577 


453.79520 


17.0510 


0.1747 


455.75470 


17.2530 


0.1547 


457.78640 


17.3390 


0.1647 


1560.39200 


15.5760 


0.0550 


1560.39600 


15.5580 


0.0550 


1561.40400 


15.5380 


0.0493 


1561.40800 


15.5300 


0.0507 


1564.40100 


15.5250 


0.0465 



Table 3 — Continued 



Event Name Pass Band time (MJD) Magnitude uncertainty 



MACHO-96-BLG-12 MACHO-Red 



MACHO-Blue 



CTIO 



MACHO-98-BLG-6 MACHO-Red 



1564.40500 


15 


,5670 


0.0409 


I joj.4UjUU 


15, 


.4540 


0.0437 


1 jOJ.4UoUU 


1 J, 


AQ^ A 
.DoDU 


A A/1 TO 

U.U4 /y 


1 jOO.Z /oUU 


1 J 


^ ^AA 


A A /I AO 


lj00.2o2uU 


1 c 

ij, 


OOA 


A A/1T1 


AAA niCiA n 


1 A 

14, 


.Z /jU 


A A1 1 O 


4j3. /j4 /U 


1 A 

14, 


.2670 


A AIOO 

U.UioZ 


4 J /. /oD4U 


1 A 
14, 


.ZooU 


A A100 


4dj.o/odU 


1 /I 
14 


.3o30 




465.65960 


14 


.2380 


0.0542 


468.75260 


14 


.3240 


0.0736 


4/i.oy /UU 


1 /I 

14 


.Zo/U 


A Aocn 
U.UZDy 


AHA n(\C\C\C\ 

4 /4. /uuyu 


1 /I 
14, 


. jdjU 




4/d.d2odU 


14, 


.2410 


0.0336 


4iSU.0l /lU 


1 /I 

14, 


.Z43U 


A A'2A1 


441. 1 yy^r\J 


ID. 


0Q'2A 

•ZooU 


A A'2Q0 


AAA '70QQf\ 

440. /ZooU 


1 c 

ID. 


.DZ3U 


A AOAA 

u.uzuu 


4D/. /iso4L) 


ID 


.Z/4U 


A AC1 1 


Alt A.nmn 
4/ l.oy /UU 


ID 


10AA 


0.0336 


476.62860 


15 


.3130 


0.0447 


480.61710 


15 


.2760 


0.0409 




15. 


.3200 


A A^ 1 1 


AQQ Aoo^n 


1 C 

I J. 


.z/zu 


A ACQA 


jUU.oUj /U 


15 


.2530 


A A/lOO 


jUl.uiZJU 


1 J 


.ZOOU 


A A'J'TO 

U.U3 / Z 


1 o34.4uUUU 


1 1 

13. 


.39zU 


A Al A/l 

u.uiy4 


10J4.4U4UU 


1 1 

Id, 




A AO 1 

U.UZIZ 


1 £'20 1AOAA 

lojy.zuzuu 


1 1 

13, 


.3890 


A AOyl'^ 

U.UZ43 


lt>4U. 14cSUU 


I J. 


.J04U 


A AOA'J 

U.UZUJ 


1640.15200 


13 


.3720 


r\ ATI'** 

0.0222 


1653.21000 


13, 


.3330 


0.0255 


1653.21400 


13, 


.3350 


0.0255 


1661.15700 


13, 


.3090 


0.0203 


1661.16100 


13, 


.3150 


0.0194 


1668.01200 


13, 


.2860 


0.0232 


1164.77190 


15, 


.8270 


0.1180 


1168.75560 


15, 


.9450 


0.1339 


1318.40870 


15, 


.8870 


0.0590 


1319.40610 


15, 


.8430 


0.0475 


1321.46010 


16, 


.1810 


0.1448 


1323.46890 


15, 


,9430 


0.1806 


1324.42410 


16 


,2700 


0.2814 


1325.41330 


16, 


.0320 


0.1448 


1325.42910 


15, 


,9010 


0.1220 


1326.46850 


16 


,2210 


0.1587 



Table 3 — Continued 



Event Name Pass Band time (MJD) Magnitude uncertainty 



MACHO-Blue 



MPS 



MACHO-99-BLG-1 MACHO-Red 



MACHO-Blue 



11/^/1 -7-7 1 nn 

i io4. / / lyu 


17. 


,1970 


0. 1856 


i loo. /JJOU 


lo. 


,9490 


0. 1966 


i jl6.4U6 /U 


17. 


, 1 240 


0.1 170 


1319.40610 


17, 


,1600 


0.0814 


1321.46010 


17, 


.1580 


0.1210 


1323.4do90 


17, 


,3430 


0.2824 




17. 


,1750 


0. 1796 


i Jzj.4zy lu 


17. 


.1360 


0. 1677 


132o.4dojU 


17, 


,2;S40 


A 1 /I AO 
0.1406 


132/.4jUjU 




,1130 


O.Ooo / 


0'2A'2 '2Q'2f\f\ 


13, 




A A'2A'2 
0.0303 


T3AC lATAA 
ZiUj.iO /{){) 


13. 


ACTA 


0.0243 


T5AO /I AAAA 

2306.40000 


13. 


,6610 


0.0243 


2318.39500 


13. 


,5660 


0.0290 


2323.41900 


13, 


,5270 


0.0266 


Tl/IO 0/;£AA 

23d2.2dd00 


13, 


,2020 


0.0437 


2430.U34UU 


13, 


OA< A 

,2050 


A AO AO 

0.0202 


2439.03100 


13, 


,2370 


A AOTO 

0.0278 


O/I /I A A10AA 

2440.03600 


13. 


on A 
,2320 


0.0243 


"> A A \ A/1 nCIA 

2441.04y00 


13. 


,2510 


0.0266 


TIAA 1 1 A1A 

230o. / 1030 


12, 


AAIA 

,9530 


A f\A <A 

0.0450 


OIAA 01 TCA 
23OD.0I 00 


13, 


Al 1 A 
,0110 


A A'2n< 

0.0395 


23OO.JJO/O 


12, 


AAOA 

,90o0 


A AAAC 

O.Oooj 


2308.65670 


13. 


.1460 


0.1479 


2308.71270 


12. 


,9020 


0.0959 


TIAO "70 1 1 A 

230o. /bilO 


Iz, 


,9300 


A Ainc 
0.0395 


T21 £ /CO 1 CA 

23iD.ooljO 




,06 /O 


0.0d3o 


oil /C *;0'2'7A 

231d.do3/0 




,91j0 


0.0535 


2328.59630 


12, 


,7660 


0.0782 


2328.77350 


12, 


,7720 


0.0593 


441.78840 


15, 


.3500 


0.0301 


442.78270 


15. 


,3320 


0.0251 


443.75950 


15. 


,3800 


0.0428 


452.77860 


15. 


,3780 


0.0504 


455.77750 


15. 


,3460 


0,0400 


459.75610 


15, 


.3650 


0.0419 


463.75080 


15, 


.3870 


0.0629 


463.75430 


15, 


.3580 


0.0600 


465.73690 


15, 


,6230 


0.1637 


466.70240 


15, 


,3400 


0.0687 


441.78840 


16. 


,2880 


0,0428 


442.78270 


16, 


.2840 


0.0327 


443.75950 


16. 


.2980 


0.0629 


452.77860 


16, 


,3700 


0.0795 
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Table 3 — Continued 



Event Name Pass Band time (MJD) Magnitude uncertainty 



455.77750 16.3130 0.0648 

459.75610 16.2700 0.0854 

463.75080 16.1810 0.0972 

463.75430 16.1740 0.1041 

465.73690 16.7720 0.3813 

466.70240 16.4000 0.1349 

CTIO 2732.31900 13.3080 0.0232 

2733.20300 13.3090 0.0266 

2769.22400 13.4460 0.0194 

2778.22500 13.5420 0.0187 

2778.99300 13.5280 0.0222 

2784.13200 13.5770 0.0222 

2789.15200 13.6390 0.0522 

2794.96300 13.7050 0.0232 

2819.99200 13.9730 0.0255 

2820.99400 14.0080 0.0232 

MPS 2688.61020 12.6050 0.0328 

2688.71570 12.5720 0.1494 

2689.47020 12.6120 0.0841 

2689.55590 12.6540 0.0465 

2689.61810 12.6300 0.0507 

2689.71450 12.8030 0.1658 

2689.82380 12.6190 0.0382 

2690.56900 12.6450 0.0493 

2690.63600 12.6540 0.0221 

2690.66750 12.6410 0.0395 

MACHO-99-BLG-8 MACHO-Red 1165.75610 11.7260 0.0194 

1166.76670 11.7530 0.0155 

1168.75920 11.7830 0.0155 

1318.41170 11.7830 0.0153 

1319.40950 11.7760 0.0153 

1323.47160 11.7490 0.0155 

1324.42880 11.7530 0.0157 

1325.41700 11.7610 0.0155 

1325.43250 11.7630 0.0155 

1326.47170 11.7480 0.0155 

MACHO-Blue 1168.75920 13.4070 0.0172 

1323.47160 13.3740 0.0177 

1324.42880 13.3880 0.0188 

1325.43250 13.3850 0.0172 

1326.47170 13.4110 0.0167 

1327.43530 13.4010 0.0157 

1329.46640 13.3660 0.0163 

1330.42570 13.3340 0.0163 

1331.46440 13.4090 0.0447 
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Table 3 — Continued 



Event Name Pass Band time (MJD) Magnitude uncertainty 



MPS 



1342.45200 


13.3370 


0.0345 


2630.39700 


11.3820 


0.0212 


2630.40000 


11.3760 


0.0194 


2630.40300 


11.3710 


0.0255 


2630.40600 


11.3960 


0.0243 


2632.39400 


11.3020 


0.0290 


2632.39600 


11.3220 


0.0355 


2632.40000 


11.3240 


0.0355 


2632.40200 


11.3840 


0.0342 


2639.37200 


11.3020 


0.0243 


2639.37500 


11.2510 


0.0203 


2688.60420 


12.4140 


0.0179 


2688.60530 


12.4340 


0.0186 


2688.60680 


12.4310 


0.0186 


2688.71320 


12.3960 


0.0493 


2689.46830 


12.4360 


0.0232 


2689.55290 


12.3980 


0.0221 


2689.61620 


12.4230 


0.0232 


2689.71320 


12.4310 


0.0202 


2690.47590 


12.3160 


0.0564 


2690.56730 


12.4250 


0.0202 



Note. — Complete data set available in electronic version. MJD = JD — 

248623.5 days. 



Table 4. Photometric Calibration Coefficients 



Event Name 


a 


b 


c 


d 


MACHO-104-C 


0.8176 


0.1824 


0.1828 


0.8172 


MACHO-96-BLG-5 


0.8076 


0.1924 


0.1804 


0.8196 


MACHO-96-BLG-12 


0.8176 


0.1824 


0.1828 


0.8172 


MACHO-98-BLG-6 


0.8191 


0.1809 


0.1829 


0.8171 


MACHO-99-BLG-1 


0.8169 


0.1831 


0.1826 


0.8174 


MACHO-99-BLG-8 


0.8188 


0.1812 


0.1829 


0.8171 
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Table 5. Microlensing Parallax Fit Parameters 



Event 


fMR 


fMB 


fcTIO 


Imps 


to (MJD) 




t (days) 


V (km/sec) 


e 






104-C 


1.00(1) 


0.99(2) 






508.3(6) 


0.15(1) 


220(2) 


77(4) 


-1.08(7) 


1.47 


1051 


96-BLG-5 


0.12(3) 


0.12(3) 


0.13(3) 




1763(1) 


0.018(6) 


2000(500) 


30.9(1.3) 


-0.84(6) 


1.58 


2395 


(HST) 


0.28(1) 


0.30(1) 


0.33 




1767(1) 


0.048(6) 


970(20) 


30.9(1.3) 


-0.87(7) 


1.59 


2371 




0.31(1) 


0.33(1) 


0.37 




1768(1) 


0.054(7) 


900(20) 


31.0(1.3) 


-0.88(8) 


1.59 


2363 


96-BLG-12 


0.87(2) 


0.89(3) 


0.90(2) 




1743.4(3) 


-0.11(2) 


294(5) 


47.5(1.3) 


-1.23(9) 


2.11 


5914 


98-BLG-6 


0.65(14) 


0.60(13) 


0.68(15) 


0.66(13) 


2388(3) 


0.16(4) 


490(50) 


79(5) 


-1.7(2) 


1.20 


802 


99-BLG-l 


0.96(9) 


0.98(10) 


1.0(1) 


0.97(7) 


2712(1) 


0.23(4) 


231(13) 


43.9(9) 


-1.85(2) 


1.54 


1706 


99-BLG-8 


0.75(12) 


0.73(12) 


0.76(12) 


0.79(13) 


2732.1(4) 


0.17(1) 


240(20) 


62(5) 


-1.53(3) 


2.34 


2280 



Note. — MJD = JD - 248623.5 days. 



Table 6. Microlensing Parallax Likelihood Mass Estimates 



Event 


location 


fMR 


t (days) 


V (km/sec) 


Confidence Levels P{M/Mq < N) 
p = 5% P = 16% P = 50% 


P = 84% 


P = 95% 


104-C 


bulge 


1.00(1) 


220(2) 


77(4) 


0.35 


0.62 


1.15 


2.2 


3.94 


96-BLG-5 


bulge 


0.12(3) 


2000(500) 


30.9(1.3) 


7.2 


12.8 


27 


69 


160 


(HST) 


bulge 


0.28(1) 


970(20) 


30.9(1.3) 


1.64 


2.93 


6.3 


15.8 


37 




bulge 


0.31(1) 


900(20) 


31.0(1.3) 


1.41 


2.53 


5.4 


13.6 


31 


96-BLG-12 


bulge 


0.87(2) 


294(5) 


47.5(1.4) 


0.33 


0.62 


1.29 


3.1 


6.7 


98-BLG-6 


bulge 


0.65(14) 


490(50) 


79(6) 


0.94 


2.6 


5.7 


12.5 


24 




SGR 


0.65(14) 


490(50) 


79(6) 


1.23 


1.61 


2.52 


4.2 


6.7 


99-BLG-l 


bulge 


0.96(9) 


219(9) 


42.9(9) 


0.14 


0.29 


0.68 


1.86 


4.6 




SGR 


0.96(9) 


219(9) 


42.9(9) 


0.10 


0.16 


0.31 


0.70 


1.51 


99-BLG-8 


bulge 


0.75(12) 


240(20) 


62(5) 


0.27 


0.56 


1.19 


2.78 


6.0 



Note. — Bold-faced type indicates the parameters that are considered to be most likely. Event 96-BLG-5 has parameters for three different fits 
listed. The first fit is the fit with no constraint on the source brightness, while the second and third fits have the lensed flux fixed to a value based 
upon our HST observations. The fit labeled HST is the best fit, while the third fit is provided to indicate the effect of the source flux uncertainty 
on the mass Umits. For events 98-BLG-6 and 99-BLG-l, mass estimates based upon bulge and SGR sources are presented. In all cases, it is most 
likely that the source star is in the bulge. 
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Table 7. Mass & Magnitude Estimates for the MACHO Microlensing Parallax Events 



Event 


M/Mq 


Mms/M© 




sep-MS 










AC/fo 


104-C 


^•^ -0.5 


0.74 


2.7kpc 


40 mas 


17.3 


3.5 


3.5 


3.5 


3.2 


96-BLG-5 




















96-BLG-12 


1 3+1-** 

-^■■'-0.7 


0.75 


2.0 kpc 


28 mas 


18.0 


2.1 


2.2 


2.2 


2.3 


98-BLG-6 


^■^-0.9 


0.88 


5.7 kpc 


5 mas 


20.1 


2.2 


1.9 


1.6 


1.1 


99-BLG-l 


7+1-2 

' -0.4 


0.40 


1.7kpc 


17 mas 


18.9 


1.8 


3.2 


3.6 


3.9 


99-BLG-8 


1 2+1'^ 


1.2 


1.6 kpc 


25 mas 


16.3 


1.3 


0.7 


-0.3 


-1.1 



Note. — These are the parameters of the "most likely" main sequence star lenses for our best microlensing parallax 
events. For MACHO-96-BLG-5, a main sequence lens is ruled out. 



